Jastrow-Correlated Wavefunctions for Flat-Band Lattices 
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The electronic band structure of many compounds, e.g., carbon-based structures, can exhibit 
essentially no dispersion. Models of electrons in flat-band lattices define non-perturbative strongly 
correlated problems by default. We construct a set of Jastrow-correlated ansatz wavefunctions 
to capture the low energy physics of interacting particles in flat bands. We test the ansatz in a 
simple Coulomb model of spinless electrons in a honeycomb ribbon. We find that the wavefunction 
accurately captures the ground state in a transition from a crystal to a uniform quantum liquid. 
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PACS numbers: 71.10.-w, 73.22.-f, 71.10.Pm 



I. INTRODUCTION 

Quantum correlated phases and rich phase diagrams 
can manifest even in systems with very little contribu- 
tion from the kinetic energy. The long-range part of the 
interaction itself can lead to interesting and unexpected 
physics. Hubbard's work on lattice models found that 
the long-range part of the Coulomb interaction alone sup- 
ports insulating crystals that stabilize at odd denomina- 
tor lattice fillings X His analysis brings up an interesting 
question: Can interactions by themselves lead to new and 
interesting physics? 

Particles hopping with an otherwise large kinetic en- 
ergy can interfere to show no dispersion in certain lat- 
tices geometries (see, e.g., Ref. 0). Such flat-bands arise 
as particles hop among a few sites and interfere to form 
localized states. At first one would expect such lattices 
to exhibit only classical localization effects. However, 
the absence of any dispersion leaves the interaction as 
the de facto dominant energy scale thus allowing inter- 
esting quantum many-body effects to appear solely from 
off-diagonal contributions in the interaction or, equiva- 
lcntly, from non-commuting lattice density operators. 

The paucity of theoretical work on interacting flat- 
band systems stems from implicit difficulties in solving 
even the simplest models. Vanishing band curvature im- 
plies very little or no screening of the interaction. Off- 
diagonal terms in the interaction can therefore be rele- 
vant in such systems. Furthermore, the absence of dis- 
persion results in a single dominant term in the simplest 
flat-band models: Hyb = Tfb^Vb' where *Pfb projects 
into a flat band and V is a two-body interaction. With 
no small parameter, conventional perturbative methods 
fail. Work on flat-band lattices has become more pressing 
with the discovery of systems hosting flat bands. 

The discovery of carbon based nanostructures, e.g., 
carbon nanotubes^ and graphene^, opened the possi- 
bility of a geometric tuning of electronic band struc- 
ture. Detailed calculations show flat bands in a vari- 
ety of physically relevant contexts, including: the edges 
of two dimensional graphene^ one dimensional graphene 
nanoribbons,— ~— hydrogenated graphene nanoribbons, 8 
collapsed carbon nanotubes,— hydrogenated nanotubes~2 
graphene dotsjii and graphene antidots 1 ^. Flat bands 



can also be found in a wealth of other compounds and 
engineered in optical lattices^ The ubiquity of flat-band 
systems imply that strong interaction effects demand fur- 
ther study. 

We propose that studies in an entirely different con- 
text, the quantum Hall regime, can shed light on the 
problem of strong correlation in flat-band lattices. Mod- 
els of the quantum Hall effects project into a single Lan- 
dau band in the high field limit. The degenerate Landau 
level leaves the Coulomb interaction as the sole term in 
the simplest models to yield direct analogies with models 
of flat-band lattices in the absence of a magnetic field. 

Parallels between quantum Hall and zero-field lattice 
formalisms have been drawn at the Hamiltonian level. 
Early wor k 14 ' 15 compared models of Laughlin's states^ 
in the cylindrical geometry^ to lattice models of spin- 
less fermions. Here it was shown that, on thin cylinders, 
components of Laughlin's wavefunction can be thought 
of in terms of Hubbard's classical crystals.— The work on 
thin cylinders marked a key advance by noting an im- 
plicit one dimensional structure in quantum Hall mod- 
els. Subsequent analyses explored intriguing aspects of 
this lattice-quantum Hall connection.— 

We argue that a connection between flat-band and 
quantum Hall formalisms can be drawn at the level of 
ansatz wave functions. Elegant but accurate many-body 
wavefunctions quantitatively model the physics of the 
quantum Hall regime . 16 i 19 i 20 Here we build a similar class 
of Jastrow-correlated wave functions^ as ansatz states 
for flat-band lattices. We numerically test a quasi-one 
dimensional example in a toy model of a flat band in 
graphene nanoribbons (Figure [T]) of up to 476 sites. We 
find that the ansatz state accurately captures the ground 
state in a remarkable transition from a Wigner crystal to 
a uniform quantum liquid driven entirely by interactions. 



II. JASTROW-CORRELATED ANSATZ 
WAVEFUNCTIONS 

We construct Jastrow-correlated first quantized wave- 
functions, applicable to flat-band lattices. We consider 
a family of states motivated by Jastrow-correlated states 
written in the lowest Landau leve h 16 ' 19 For a generic two- 
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body interaction that depends only on relative separa- 
tion, i?FB commutes with the center of mass operator. 
We may therefore construct states of the form, V&cM^relj 
where fcM and are functions of the center of 

mass and relative coordinates, respectively. We propose 
that the following set of unnormalized relative-coordinate 
wavefunctions are energetically favorable ansatz states 
for a broad class of flat-band lattice models: 

l*rel>= II (4a-4,a) TO V">. (1) 
j<k,{A} 

where the symbol {A} denotes a set of one dimensional 
(ID) chains that cover a bravais lattice. Each bond- 
oriented chain uniquely covers bonds and the two end 
sites lie at the lattice edges. The first quantized operator 
for particle j, (a| A ) n , creates a state localized at the po- 
sition r ^ = 5 \ + nh\ , where n € N and 5 a is the position 
vector for a unique edge starting site of the chain A (e.g., 
S = for a ID lattice and {6} = {(n° x , n° y ), (n° x ■ n° y = 0)} 
for a square lattice). The set of bond vectors bA point 
from the starting site to a nearest neighbor site along the 
chain. (This chain covering procedure must be modified 
for some non-bravais lattices.) We define real space single 
particle basis states as (r|(a A )™|0 >= \/n!/ n (7)w(r— r^), 
where w are Wannier functions. / is a variational func- 
tion. For the model studied here we choose f n (j) — 
exp(7 2 n 2 /2)/-\/nT, where 7 is a variational parameter. 

The Jastrow factor in Eq. (flj enforces Gutzwiller pro- 
jection. It attaches an integer number of correlation 
holes, m c , to each particle in the wave function ip" at a 
filling, u*, to thereby form a state at the reduced filling 
v = v* /(m c v* + 1), akin to a procedure introduced in 
the quantum Hall regimefiS Here v refers to the number 
of particles per basis state. In the following we consider 
fcrmions. The choice m c = 2p with antisymmetric ip v 
preserves the antisymmetry of Viol- 
in the following we test a ID example of Eq. ([1]) and 
a specific choice for if)" . We consider a Hartree-Fock 
state, r* =1 = Det[(a])™]|0) = \[ j<k {^ - 4)|0>. With 
this choice the relative coordinate wavefunction takes a 
form similar to the Laughlin state : 16 ' 22 

/ , , s 2p+l 

i^D- re i)=n( a ]-4) io>. (2) 

j<k 

For an A^-particle system, the number of particles per 
basis state is then v = N/[(2p + 1)N — 2p\. We focus on 
p = 1, which corresponds to v = 1/3 for N — > 00. 

The variational parameter, 7, tunes the above wave- 
functions between distinct limits. For 7 — > 00 the wave- 
functions describe crystals (Wigner crystals) because the 
basis states are highly localized. For intermediate val- 
ues, 7 ~ 1, we have a charge density wave. For 7 — > 
we obtain uniform liquid states. To see this one notes 
that the Wannier functions can be written in terms of 
polynomials in Fourier transform spaced For a ID lat- 
tice along x we have w(x — x n ) — J- exp(iqn)w(q) , where 



J- denotes the Fourier integral over all q-space and w(q) 
is the Fourier transform of the Wannier function at the 
edge. Then Eq. [5] becomes a Laughlin-like state: 

iini*r D . rcl =^n( ei93 " em ) 2p+1 n*fe)' ( 3 ) 

j<k j 

which describes a uniform quantum liquid. 

III. MODEL OF INTERACTING FERMIONS IN 
A FLAT BAND 

We test the validity of Eq. @ as an ansatz ground state 
for a model of N spinless fermions interacting through a 
truncated Coulomb interaction in a flat band of a honey- 
comb ribbon. We consider a second quantized flat-band 
model: 

where the matrix elements Vt n \ are determined by the 
form of the interaction V(|r — r'Q and flat-band basis 
states (j) as: 

V {n} = J drdr'<(r)^(r')V(|r-r / |)^3(r')^uW.(5) 

The second quantized operator creates a fcrmion in 
the state <fi n . 

Flat-band basis states follow from the lattice structure. 
In the following we work in units of the honeycomb bond 
length, a = 1, and Coulomb energy, e 2 /eao (~ 10 5 K 
for a carbon bond in vacuum). We consider a quasi- 
1D honeycomb ribbon of width L y = y/3N y along the y 
direction. Along the ribbon we allow the range of x to be 
infinite, but restrict the centers of local electron orbitals 
to lie in N x unit cells along a length of L x = 3N X . For 
simplicity, we assume Gaussians of width ~ /? localized at 
each site of the honeycomb lattice. The hopping energy 
between sites can be used to determine f3 (or vice versa) 
prior to projection into the flat-band. In the honeycomb 
geometry the flat-band basis states are: 

2Ny 

0„(r) =i^]Ti 5 . e -^ R ") 2 /^ 2 , (n = l,..,N x ), (6) 

where Np is a normalization constant, = (X*, Y s ) is 
the location of the sth node in the nth unit cell, and the 
node factor t s has the form of t2 U -i = <2« = (— 1) " +1 for 
u = 1, .., N y & The node locations arising from the under- 
lying honeycomb structure are: R* = (3n — e s , \fi{u — 
1/2)), where e s = 2 (e s = 1) for s = 2u - 1 (s = 2u) 
and u = 1, ..,N y . In the following we set N y = 3. Us- 
ing the above Wannier state 4> n as a single-particle basis 
state localized at the nth unit cell, we can construct the 
interaction matrix elements analytically. 



We stress that the model considered here idealizes flat- 
bands of armchair nanoribbons. The Gaussian localized 
states chosen in Eq. |6] are only convenient approxima- 
tions for orbitals of electrons in flat bands of suspended 
graphene nanoribbons. For example, the 7r z Slater or- 
bitals of carbon atoms could be used instead.- These or- 
bitals are, in contrast, exponentially localized and loosely 
correspond to /3 ~ 1. We expect less localization (larger 
ft) in carbon structures with adsorbates (e.g., hydrogen^). 
It is straightforward to replace the Gaussian orbitals in 
our formalism with density functional theory results for 
specific systems to make concrete predictions for experi- 
ments. Furthermore, we have ignored effects from nearby 
bands. 

We now compute the interaction matrix elements 
in the flat-band basis with Gaussian orbitals. The 
Fourier expansion of the interaction is V(r) = 
e'yQ J dQV{Q)e^ x /TT with Q = tt/Vs and V(Q) = 
( e2 / e ) Em VmL m (Q 2 +Qo), were the L m are the Laguerre 
polynomials, m is an odd integer for spinless fermions, 
and V m = v / 7r(2m — l)!!/(2 m+1 m!) is a pseudopotential 
parameter for the full Coulomb interaction.™ We choose 
only the m = 1 term, which corresponds to an effective 
short-range interaction V(r) ~ V 2 <5(r). The matrix ele- 
ments are then: 

{«} 

where we define: 

T {s} = ViT aM T^J(VtoN$), 

t ss , = t s t s , [ Lv d ye M°y-((y- YS ) 2 +(y- Ya 'fyw*)\ 

Jo 

r nin4 = (X£-X£)/(20), 
r n2 n 3 = (X£-X%)/(20), 
AR = (X°\ + X£ - X£ - X%)/(20), and 

" = P~ 2 + Ql (8) 
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FIG. 1: (Color online) Greyscale plot of numerically cal- 
culated solutions for the electron density of the flat-band 
Coulomb problem in a honeycomb ribbon with N — 12 elec- 
trons and N x ~ 34 unit cells zoomed into 12 central unit 
cells. Sites sit at vertices in the depiction of the underlying 
lattice. The three graphs show three distinct states that arise 
from increasing width parameter of basis states: (a) Wigner 
crystal (/3 = 2.0), (b) Charge density wave (/? = 3.5), and (c) 
Uniform state (j3 = 5.0). 
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FIG. 2: (Color online) The local density of each unit cell, p n , 
plotted versus the unit cell index along the honeycomb ribbon 
with N = 12 particles. The circles, diamonds, triangles, stars, 
and squares are obtained from width parameters /3 = 2, 3, 4, 5, 
and 6, respectively. The lines are a guide to the eye. 



IV. NUMERICAL TEST OF ANSATZ 

We test our ansatz by solving Eqs. with numer- 

ical diagonalization on finite-sized ribbons. We restrict 
our ground state to the Hilbert space with the center of 
mass at the center of the ribbon. In the thermodynamic 
limit the center of mass is a good quantum number but 
in finite sized systems we omit edge effects in extreme 
unit cells by restricting our basis. We test several dif- 
ferent system sizes and show below that our results have 
converged with increasing system size. 

Figure []] plots a zoom in of the quasi-lD charge density 
over a few central unit cells of an N = 12 and N x = 34 
system computed from our fiat-band model, Eqs. ((U)-©. 
The charge density is p(r) = J2n=i \4 > n{ r )\ 2 Pnt where 
p n =< ^fbIc^CjjI^fb > is the local density of the nth 



unit cell and ^fb is the exact ground state. The density 
is plotted against the underlying honeycomb lattice to 
show that each unit cell consists of 14 sites in this exam- 
ple. Figure [1] shows the evolution from a Wigner crystal 
configuration (a), to a charge density wave (b), and then 
to a uniform state (c) as we increase the extent of basis 
states. The same features are exhibited in Figure [21 the 
local density plot for each unit cell. 

We check the validity of Eq. (J5J) in describing the 
ground state of Eqs. (HJ-® by taking overlaps. We gen- 
erate Eq. @ numerically by noting that it is the exact 
eigenstate of a model of the lowest Landau level of elec- 
trons on a cylinder ! 15 ' We obtain the maximum overlap 
with the choice 7 = 3//3. The right axis of Figure [3] plots 
the overlap of Eq. © with the ground state of Eqs. flU)- 
©. The left axis plots the deviation in the local den- 



4 




FIG. 3: (Color online) The left axis plots the deviation in local 
density between neighboring unit cells at the sample center 
versus the basis state width parameter for three different sys- 
tem sizes . The right axis plots the overlap between the ansatz 
state, Eq. ([2|, and the numerically computed ground state of 
Eqs. dm-©. 
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FIG. 4: (Color online) The excitation energy plotted versus 
the basis state width parameter for three different system 
sizes. Open (closed) symbols indicate the energy gap in the 
relative (center of mass) degrees of freedom. 

sity from the center to the nearest neighboring unit cell, 
Ap c = \p(N x /2) ~ P(N x /2-i)\- A/9 C = 1 and correspond 
for a crystal and a uniform state, respectively. We find 
that the overlap is essentially constant and larger than 
99.61% for all (3 thus showing that the ansatz wavefunc- 
tion accurately captures the ground state of the trun- 
cated Coulomb interaction in all parameter regimes. We 
expect the overlap to decrease as we increase the range of 
the interaction or change the form of the Wannier func- 
tion. Overlaps can be improved by choosing a better 
variational form for f('j). 

It is interesting to note that the ansatz state adiabati- 
cally connects the Wigner crystal to a uniform, Jastrow- 
correlated liquid state based on the Laughlin wavefunc- 
tion. This crossover was pointed out in the quantum Hall 
context in Rcf. 15. Here we have shown that a mapping 
of the quantum Hall state to a zero- field fiat-band lattice 
preserves this crossover. 

Figure @] plots the energy gap as a function of /3. In 
an infinitely large system, excitations can either leave 



the center of mass intact (excitations in the relative co- 
ordinate degrees of freedom) or shift the center of mass 
(center of mass excitations). We find a small but finite 
gap for the uniform state (large /?) in the relative degrees 
of freedom. The gap appears to remain robust against 
system size suggesting that the gap we have computed is 
indeed a bulk gap, thus signaling the intriguing possibil- 
ity of an incompressible uniform state. We note, however, 
that the center of mass gap becomes comparable to the 
relative coordinate gap for large /3. It is possible that ex- 
citations in the relative degrees of freedom are spatially 
large in this case and therefore sensitive to edge effects. 
Periodic systems will allow a more accurate estimate of 
the bulk gap. 



V. SUMMARY 

We have constructed a general class of Jastrow- 
correlated wavefunctions applicable to flat-band lattices. 
We studied a toy model of spinless fermions in a fiat- 
band in graphene nanoribbons in the absence of a mag- 
netic field. We find that the wavefunction accurately 
captures an intriguing transition from a Wigner crystal 
to a Jastrow correlated liquid, driven entirely by inter- 
actions. Our method can be generalized to accurately 
model graphene nanoribbons and other flat-band lat- 
tice structures^— and can be adapted to study mag- 
netism recently explored in experiments on graphene 
nanoribbons.— 

The Jastrow-correlated wavefunctions constructed 
here are versatile. They can be combined with density 
functional theory estimates of Wannier functions to si- 
multaneously capture band structure and the effects of 
strong interactions. They can be generalized to apply 
to a wide variety of lattices, describe correlated bosons, 
or incorporate spin. For example ^ TC \ can quantify a 
connection between Gutzwiller-RVB and quantum Hall 
states^ Gutzwiller projected d-wave spin-singlet paired 
states in a two dimensional square lattice can be con- 
structed by writing i\) v as an antisymmetric product of 
singlet paired states, oc -4[g(i"i — r 2 )g(r 3 — r4)...]. Here 
an effective pairing among i\) v * particles could arise from 
an over-screening of long range repulsion, as in the quan- 
tum Hall regime) 20 i 27 ' 28 but in flat or narrow-band lattice 
models. 
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